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ABSTRACT 



Aims. We aim to characterize the (/-band variability of young brown dwarfs in the Taurus Molecular Cloud and discuss its origin. 
Methods. We used the XMM-Newton Extended Survey of the Taurus Molecular Cloud, where a sample of 1 1 young bona fide brown 
dwarfs (spectral type later than M6) were observed simultaneously in X-rays with XMM-Newton and in the t/-band with the 
XMM-Newton Optical/UV Monitor (OM). 

Results. We obtained upper limits to the (/-band emission of 10 brown dwarfs ((/> 19.6-20.6 mag), whereas 
2MASS J04141 188+281 1535 was detected in the (/-band. Remarkably, the magnitude of this brown dwarf increased regu- 
larly from (/~19.5 mag at the beginning of the observation, peaked 6 h later at U~ 18.4 mag, and then decreased to (7 ~ 1 8.65 mag in 
the next 2h. The first OM (/-band measurement is consistent with the quiescent level observed about one year later thanks to ground 
follow-up observations. This brown dwarf was not detected in X-rays by XMM-Newton during the OM observation. 
Conclusions. We discuss the possible sources of (/-band variability for this young brown dwarf, namely a magnetic flare, non-steady 
accretion onto the substellar surface, and rotational modulation of a hot spot. We conclude that this event is related to accretion from 
a circumsubstellar disk, where the mass accretion rate was about a factor of 3 higher than during the quiescent level. 
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1. Introduction 

Young brown dwarfs (BDs) share similar properties with 
pre-main-sequence low-mass stars, suggesting that they 
are also undergoi ng a T Tauri p hase. BDs have accretion 
disks (e.g., 



Mohantv et al. 



Javawardhana et all I2003I 



2005) 
well 



and photometric 
2004a| b|). as well as spectroscopi c variability ( Schol z et al 
2005; Scholz & Javawardhana 2006), pointing to variable mass 
accretion. Besides, BD s also have an active corona, em itting 
soft X-rays (e.g.. lNeuhauser et al.ll999tlPreibisch et all2005l) . 

The gas accreting from the circumstellar disk onto the stel- 
lar surface of T Ta uri stars produces a UV/blue excess emis- 
sion (see review by Bouvier et al. 2006J). The (/-band luminos- 
ity of T Tauri stars is d irectly related to the accretion luminosity 
<GullhringetaljH99l . However, due to their faintness in the 
optical, BDs are poorly studied at short wavelengths ( Gizi s et aTl 
l2005h . During the XMM- Newton Extended Survey of the Taurus 
Molecular Cloud (XEST; Giidel et al. 2006), the X-ray emission 
of a sample of 16 bona ride BDs (spectral type later t han M6) 
was surveyed with XMM-Newton dGrosso et alJl2006bl) . Eleven 
of these BDs were located less than ~8.5' — 12' away from the 
prime target of 9 XMM-Newton pointings, which allowed for 
simultaneous observations with the XMM-Newton Optical/UV 
Monitor (OM; Mas on et alJ2001l) . that we report here. 

We present the constraints on the (/-band photometry of 
these BDs obtained with the OM in Sect. [5] In Sect. [2 we inves- 



lNattaetal.1 f2 004; 
(Scholz &Eisl6ffe 
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tigate the variability of 2MASS J04141 188+281 1535 (shortened 
to 2MASS J0414 henceforth) in the (/-band and in X-rays. We 
discuss the origin of the U -band event observed from this BD in 
Sect.EI 



2. The OM observations of BDs in the Taurus 
Molecular Cloud 

Table Ogives the list of the 11 BDs of the Taurus Molecular 
Cloud (TMC) surveyed with the OM in 9 XMM-Newton point- 
ings. About half of these XMM-Newton pointings - namely 
XEST-02, 05, 11, 14 - were made with the OM in full-frame 
imaging mode, which provides an image of the whole OM field 
of view (17'xl7') with 2x2 binned (low) spatial resolution 
(1" x 1"). The other XMM-Newton pointings - namely XEST- 
01, 04, 08, 15, 20 - were made with the OM in the default 
imaging mode, which consists of a sequence of 5 exposures 
with 5 imaging windows - allowing us to cover 92% of the to- 
tal (17'xl7') OM field-of-view - plus a small central imaging 
window (1.7'xl.7')- In each of the 5 exposures, one of the 5 
large imaging windows covers a large fraction of the field-of- 
view with 1" x 1" spatial resolution, while the small imaging 
window ensures continuous monitoring of the prime target at 
the center of the field-of-view with 075 x 075 spatial resolution. 
We indicate in Col. (7) of Table^the duration of one single full- 
frame/imaging window exposure and in Col. (8) the total number 
of similar exposures obtained on the BD. 
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Table 1. The sample of the 11 young BDs, w ith spectral type later than M6, surveyed with the XMM-Newton Optical/UV 
Monitor (OM) in the XEST l lGudel et ai] l2006). Columns (1) and (2) give the BP names and the 2 MASS c ounterparts respec- 
tively. Column (4) is the refere nces of the spect r al typ es in Col. (3), with B02=Briceno et al. (2003), G06= lGuieu et alJ (12 006 ). 
L06= Luhmanl ( I2006I) . an d M05=lMuzerolle et all ( 120051) . Column (5) gives the disk accretion rate when available ( Mohanty et al. 
120051 for 2MASS J0414.lMuzerolle et all2005l) . Column (6) is the XEST field number, plus the source number if there is an X-ray 
detection ( iGrosso et alJl2006bl) . Columns (7) and (8) indicate the duration of one single full-frame/imaging window exposure and 
the total number of similar exposure obtained on the BD. Column (9) gives the U -band magnitude or the U -band limiting magnitude 
of each window exposure for a 3<x point-source detection (this work). For 2MASS J0422, the limiting magnitude was computed for 
the broader UVWl filter (200^400 nm), which includes the f/-band. 



Name 2MASS SpTyp Ref. logM XEST Window exposure N U 

M Q yr~' ks mag 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


2MASSJ0414 


J04141 188+281 1535 


M6.25 


M05 


-10 


20 


2.5 


10 


19.5 


2MASSJ0421 


J042 15450+26523 15 


M8.5 


L06 




11 


5.0 


5.0 


>20.4 


2MASSJ0422 


J04221332+1934392 


M8 


L06 




01-062 


3.5 


4.0 


>19.9 


KPNO-Tau4 


J04272799+2612052 


M9.5 


G06 


-11.1 


02 


5.0 


4.6 


>20.4 


KPNO-Tau5 


J04294568+2630468 


M7.5 


B02 


n.a. 


15-044 


1.3 


3.9 


>19.8 


KPNO-Tau6 


J04300724+2608207 


M9 


G06 


-10.8 


14 


5.0 


2.8 


>20.6 


KPNO-Tau7 


J043057 18+2556394 


M8.25 


B02 


-11 


14 


5.0 


2.8 


>20.6 


CFHT-Tau 5 


J04325026+2422115 


M7.5 


G06 


n.a. 


04-003 


1.5 


4 


>19.8 


KPNO-Tau9 


J04355143+2249119 


M8.5 


B02 


n.a. 


08 


1.6 


4 


>19.6 


CFHT-BD-Tau 2 


J04361038+2259560 


M7.5 


B02 


n.a. 


08 


1.3 


8 


>19.7 


CFHT-Tau 6 


J04390396+2544264 


M7.25 


G06 




05-005 


5.0 


3.7 


>20.5 



In a nutshell the OM detector is a micro-channel plate 
(MCP) intensified CCD detector operated as a photon-counting 
instrument. First, an optical/UV photon collected by the 30 cm- 
diameter telescope mirrors hits a photo-cathode located at the 
backside of the detector entrance window. The photoelectron is 
then amplified by two MCPs. The resulting electron avalanche is 
converted to a photon splash by a fluorescent phosphorus layer 
at the top of a fiber taper, which illuminates the CCD pixels. 
Finally, the onboard software determines the accurate detector 
position of the photon splash's centers in each readout frame us- 
ing a centroiding algorithm, which subsamples each CCD pixels 
by 8 x 8 in-memory pixels (about 0'.'5 x 0f.'5 on the sky). 

We use the PERL script omi chain of the XMM-Newton 
Science Analysis System (SAS; version 6.5) with default parame- 
ters to run the OM imaging mode pipeline. We corrected OM sky 
source positions in the observation source list (*OM*SWSRLI* 
fits files) of residual (a few arc seconds) boresight errors by us- 
ing the 2MASS positions of OM sources as references. Only 
one BD was detected with the OM in the t/-band (305-390 nm, 
■leff = 344 nm). We made a visual inspection of the stacked sky 
image covering the full OM field produced by the pipeline to 
look for possible faint BD missed by the OM source detection 
algorithm in individual imaging window exposure, but we found 
no other detection. 

For each non-detected BD, we compute, using an IDL pro- 
cedure, the limiting magnitude of each window exposure for a 
3cr point-source detection at the detector position of the BD be- 
cause this local value is not provided by any SAS command. We 
first need to determine the relation between sky coordinate and 
detector coordinate. We rotate the fits sky image of the exposure 
window (*OM*SIMAG* fits file) using the roll angle of space- 
craft pointing, we use details of window configuration (MODES 
extension of *OM*IMAG_* fits file) to compute absolute po- 
sition within OM frame, and we apply OM distortions (intro- 
duced mainly by the fiber taper) 1 . Comparison with the obser- 
vation source list, which provides both sky and detector source 
coordinates, helps to fix small systematic offsets in our coor- 



1 Provided by the POLYNOMJV1AP2 extension of the current cali- 
bration file OM_ASTROMET_001 l.CCF. 



dinate transformation. We apply it to BD sky coordinates with 
the boresight error added to find the proper BD detector coordi- 
nates. Inside a 6"-radius disk centered at this position we esti- 
mate the sky count rate on the detector using a 2<x-clipping pro- 
cedure 2 . We check the absence of bad pixels in the aperture using 
the QUALITY extension of this fits file. We compute the corre- 
sponding count rate error at the 3<x level assuming a Poisson 
noise. We derive the count rate error of incident photons taking 
into account coincidence-loss, deadtime, and temporal sensitiv- 
ity degradation correction (see Appendix lAl. Finally using the 
OM zero point in the [/-band filter (ZPTU - 18.26 mag), we 
convert this incident count rate to U -band magnitude. Table [2 
gives the upper limits in the {/-band. We conclude that the 10 
non-detected BDs have C/>19.6-20.6mag. 

This sample of 10 non-detected BDs is comp osed nearly 
equal ly of accreting and non-accreting sources (see Grosso et al. 
2006b). The former have f/>20.4mag, and an average visible 
extinction of 1 .5 mag, which leads, after correction of the ex- 
tinction ( Carde|hetailT989l with/? v =3.1), to {/>19.4mag and 
the corresponding [/-band luminosity 3 , L\j S5 5.4 x 10 29 ergs _1 . 
Assuming that the f/-band luminosity is directly related to 
the accretion luminosity as observed for the T Tauri stars 
(Gullbrin g et al.lll998l) . we can derive an estimate of the upper 
limit of the disk accretion rate. Taking typical values for the BD 
mass (~0.07 M ), substellar radius (/?* ~ 0.5 R Q ), and disk inner 
radius (-4.5 R+) leads to M ^ 1 x lO" 1() M yr _1 , which is con- 
sistent with the disk accretion rates obtained from spectroscopy 
for some of these BDs (Table [Q. We focus on the only BD de- 
tected in the {/-band in the following sections. 



3. The variability of 2MASS J041 41 1 88+281 1 535 

Figure[T]shows the OM image of 2MASS J0414, the young BD 
detected in the t/-band with the OM. This BD is located only 24" 

2 See IDL DAOPHOT-type photometry procedures: aper.pro and 
meanclip.pro available at http://idlastro.gsfc.nasa.gov . 

3 We adopt a distance of 140 pc for the TMC and use the conversion 
factor from the OM count rates to fluxes for M-type stars available at 

http : / /xitim . vilspa . esa . es/sas/documentation/watchout/uvf lux . shtml 
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Fig. 1. The XMM-Newton Optical/UV Monitor (OM) image in 
the {/-band of the young brown dwarf (BD) 2MASS J0414. This 
sky image is the sum of ten OM central window exposures; the 
total on-source exposure was about 7 hours. The image size is 
2' x 2' with North up and East left, with pixels of 1" x 1"; the 
stretch of the color scale is logarithmic. This image is centered 
on the weak-line T Tauri star V773 Tau located only 24" away 
from the BD. The Classical T Tauri star FMTau is also visi- 
ble. The 6"-radius circle centered on 2MASSJ0414 indicates 
the aperture used in the photometry measurement. The FWHM 
of the OM in the {/-band is -176. 



away from the weak-line T Tauri star V773 Tau, the prime tar- 
get of this XMM-Newton observation located on-axis 4 . Internal 
reflection of the light of V773 Tau (£/ =13.1 mag) within the de- 
tector entrance window produced an out-of-focus ghost image, a 
"smoke ring" centered -40" South- West away from V773 Tau. 
The default 7"— 12.5" radii annulus centered on the source used 
in omi chain to estimate the sky photometry overlaps the smoke 
ring. Therefore, the estimate of the sky level is overestimated, 
and hence the photometry of 2MASS J0414 provided by this SAS 
task is underestimated. The 6" -radius disk aperture used to esti- 
mate the source+sky photometry also partially overlaps the bor- 
der of the smoke ring. The default aperture photometry needs to 
be refined. 5 

We take the sky region on the detector image (here the small 
central imaging window; *IMAG_0000*) as close as possible to 
the source. Indeed, the FWHM is about 176 in the {/-band, and 
the PSF wings are limited for this weak source, consequently 
there is no obvious signal from the source beyond 4". Therefore, 
we tune the sky annulus from4"-6" pixels, i.e., as close as pos- 
sible to the source center. This operation is not allowed in the 
SAS command omsource, and we had to use our own IDL pro- 
cedure for this aperture photometry. We keep the canonical 6"- 
radius for the source+sky aperture because this is the size that 
was used in the calibration of the coincidence-loss. Moreover, 



4 A limitation in the onboard pixel subsampling algorithm introduces 
a modulo 8 pattern, which is removed by the data reduction pipeline, 
but residual modulo 8 background fluctuations are still visible around 
this bright source. 

5 The OM light curve of 2MASS J0414 presented previously in con- 
ference proceedings l Grosso et al. 2006a) was obtained with SAS ver- 
sion 6. 1 and default aperture photometry, and was then strongly affected 
by the smoke ring. 
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Fig. 2. The OM light curve of 2MASS J0414 in the t/-band. The 
image sequence shows enlargements (field-of-view of 12" x 12" 
with North up and East left, and with 0.5" x 0.5" pixel-size) 
of the OM central small window exposures (2480 s) in the U- 
band; the stretch of the color scale is linear. The OM observed for 
-6 hours showed an increase in the brightness (At/ 1 mag). 

Table 2. Optical CCD photometry of 2MASSJ0414 obtained 
with the 1.3-m McGraw-Hill telescope at the MDM observatory 
on October 25, 2005. Filter values are the average of two obser- 
vations obtained during the same night (05 h 19 m -05 h 34 m UT and 
ll h 46 m -ll h 55 m UT). 



Name 


U 


B 


V 




mag 


mag 


mag 


2MASSJ0414 


19.49±0.25 


19.13+0.10 


19.56+0.07 



taking a smaller radius would imply the use of the PSF calibra- 
tion (note that the PSF shape is flux dependent due to the non- 
linearity of the detector), and hence the introduction of other cal- 
ibration uncertainties. We measure the count rates of sky+source 
and sky apertures on the detector, we correct the sky count rate 
from the difference in aperture area, and we apply non-linearity 
corrections (AppendixlAl to both detector count rates to retrieve 
the count rate of sky+source and source incident photons. The 
count rate of source incident photons is simply the difference 
between these two numbers. 

The resulting photometry of 2MASSJ0414 is shown in 
Fig. |2 The location of this BD on the small central imaging 
window of the OM provided a nearly continuous coverage dur- 
ing -8 h with 42 min-exposures. The magnitude of the BD in- 
creased regularly from f/~19.5 mag at the beginning of the ob- 
servation, peaked 6 h later at U~ 18.4 mag, and then decreased to 
U~ 18.65 mag in the next 2 h. 

We obtained, about one year after the OM observation, 
a follow-up observation of 2MASSJ0414 with the 1.3-m 
McGraw-Hill telescope at the MDM observatory. This source 
was observed twice, 6.4 h apart on October 25, 2005. The source 
had a constant magnitude of U — 19.49 ± 0.25 mag (Table[2}- 
The first OM U -band measurement is fully consistent with this 
ground photometry; therefore we consider this U -band level as 
the quiescent level of 2MASS J0414. 
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Fig. 3. X-ray image of V773 Tau in the 0.5-2 keV energy band 
obtained with XMM-Newton/EPIC (pn+MOSl+MOS2). The 
field-of-view and orientation are the same as that in Fig.^ The 
stretch of the color scale is logarithmic. 2MASSJ0414 is lo- 
cated on the PSF wings of V773 Tau, where no X-ray coun- 
terpart was found with the source detection algorithm. The 6"- 
radius circle indicates the region used for the selection of the 
background+source X-ray events. The 12"-width annulus cen- 
tered on V773Tau indicates the region used for the selection 
of the background X-ray events; the circular region centered on 
FMTau, another X-ray source, was excluded from this annulus. 
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Fig. 4. XMM-Newton/EPIC pn X-ray light curve of 
2MASSJ0414. The solid lines show the source+background 
light curve at the location of 2MASSJ0414 (see the definition 
of selected regions in Fig.©. The dashed/dotted lines show the 
estimate of the background level inside the extraction region for 
comparison. The time bin size is one hour. The scale of the time 
axis is identical to the one used in Fig. [2] 2MASS J0414 is not 
detected in X-rays. 



Figure [3] shows the simultaneous XMM-Newton/EPIC 
(pn+MOSl+MOS2) image of V773Tau. 2MASS J0414 is lo- 
cated on the highly structured PSF wings of this bright X-ray 
source. No X-ray source was detected at this location with the 
SAS source detection algorithm. The upper limit on the X-ray lu- 
minosity in the 0.5-8 keV e nergy range is L x < 1 .6xl0 28 erg s" 1 , 
and log(L x /iboi) ^ -3.6 (Gros soet alJl2006bl) . For compari- 
son, BDs detected in X-rays in the TMC have log(Lx/£boi) ~ 



-4.5 to -3.0 (Neuhaus CTetai]IT999t iMokler & Stelze J 120021 
iGrosso et alJl2006bT) : therefore this upper hmit is not so strong. 
We also look here for possible X-ray variability, which could 
help to detect this BD during a limited time interval. Figure 0] 
shows the source+background X-ray light curve for XMM- 
Newton/EPIC pn at the position of 2MASS J0414, compared 
with the estimate of the background level inside the 6"-radius 
extraction region. No X-ray excess is found above the back- 
ground, hence 2MASS J0414 is definitively not detected in X- 
rays. 



4. Discussion 

2MASS J0414 was identified as a (young) BD of the TMC by 
iLuhmanl (12001 . It has a visual extinction of 1.1 mag; an M6.25 
spectral type, which corresponds to an effective te mperature 
of 29 60 K; and a bolometric luminosity of 0.015 L G jLuhman 

2004) . From these values the BD radius can be estimated t o 
0.5 R . Using the evolutionary tracks from Bar affe et alJdl998l) . 
the BD mass is 0.07 M . The only published optical spectrum 
of 2MASS J0414, obtained on December 15, 2003, exhibits a 
strong excess of continuum at short wavelengths plus a strong 
Ha emission (EW(Ha)=250A: iMuzerolle etalJ 120051) in dicat- 
ing accretion from a circumsubstellar disk (Luhman 2004), with 
an es timated accretion rate of ~10 _10 M G yr _1 (Muzerolle et alJ 

2005) . Blueshifted absorption components superimposed on the 
Ha accretion emissio n line provide evidence for an outflow 
( Muz erolle et alJ2005l) . 

2MASSJ0414 has B — V 0.4 (Table©, which is ex- 
tremely blue for an object with an M6.25 spectral type. This 
strong blue excess indicates that this BD is actively accret- 
ing. Following the method used in Sect. |2 we find that its 
quiescent level, U = 19.5 + 0.2 mag, corresponds to M = 
1. 9 + 0.4 x 10~ 10 M o yr~ 1 , which is consistent with the value 
of M uzerolle et al.l J2005I) derived from spectroscopy. Therefore, 
the high disk accretion rate of 2MASS J0414 is likely the reason 
why we detected only this BD in the U -band. 

We consider 3 possible sources of U -band variability for this 
young BD, namely: a magnetic flare, non-steady accretion onto 
the substellar surface, and rotational modulation by a hot spot. 
We discuss these 3 possibilities in turn. 

XMM-Newton observations of coronal flares from active 
stars have shown an extremely good tem poral corre lation be- 
tween UV and X-ray variations ( Mitra-K raev et al] 120051) . A 
similar result was found in the optical {/-band and X-rays 
in an XMM-Newton observation of Proxima Cen, where sev- 
eral X-ray weak flares were preceded by an optical burst 
( Giid el et alJ |20 04). We can predict the flare luminosity in- 
crease in the X-rays with the UV-X-ray relationship in flares of 
Mitr a-Kraev et alJ ( 120051) . The increase from quiescent to flare 
peak in the [/-band is equivalent to a spectral luminosity den- 
sity of 2.4 x 10 26 erg s" 1 A" 1 . The Mitr a-Kraev et alJd2005h rela- 
tionship, based on the UV Wl OM filter (lEhle et alJ2005l) . yields 
~ 1 x 10 29 erg s _1 in the 0.5-8 keV energy range, which is a lower 
limit on the X-ray luminosity increase because the U -band filter 
covers roughly half of the wavelength range of the UVWl fil- 
ter (200^1-00 nm). Therefore, a UV flare would produce an X- 
ray flare at least ~10 times brighter than the upper limit on the 
X-ray luminosity of 2MASS J0414, hence it would be easily de- 
tectable. Yet, we do not detect 2MASS J0414 in X-rays during 
the U -band maximum of luminosity. Moreover, the luminosity 
increase lasts for at least 6 h and exhibits a slow continuous rise. 
In contrast, the f/-band flares observed from Proxima Cen are 



N. Grosso et al.: A U-band survey of BDs in the TMC with the OM 



5 



impulsive , and the duration of the longest flare is shorter than 
one hour ((iiidel et alJ20 04). We thus conclude that the U-band 
luminosity rise is unlikely to result from a magnetic flare. 

Emission line variability occurrin g on timescales of hours 
has be en reported for young BDs by IScholz & Javaward hana 
(2006). The line variability is interpreted as non-steady accre- 
tion onto the substellar surface with the accretion rate varying 
by one order of magnitude over a few hours. The {/-band lu- 
minosity rise observed for 2MASS J0414 could similarly result 
from an accretion burst onto the central object. 

Assuming a purely radial free-fall motion of the gas from 
the disk's truncation radius towards the corotating BD, the 6 h 
timescale of the f/-band rise would correspond to the travel time 
from a distance of 3.5 substellar radii (i.e., 1.8 R Q ) above the sub- 
stellar surface 6 ; the terminal velocity would be about 200 km s _1 . 
The Keplerian period at the disk's truncation radius is then about 
35 h, whic h is in the range of rotational periods observe d for BDs 
(4-1 lOh Scholz & Eisloffel 2005; Zapatero Osorioeta0[2006). 
The U -band increase of ~ 1 . 1 mag would then translate into a 
mass accretion rate enhanced by a factor of nearly 3 compared 
to the quiescent level, assuming that the [/-band luminosity is 
directly relate d to the accretion lu minosity as observed for the 
T Tauri stars dGullbring et alll998l) . 

As an alternative to the accretion burst, the 2MASSJ0414 
light curve may simply reflect the rotational modulation of the 
object's luminosity by a hot accretion shock located on the sub- 
stellar surface. Evidence for rotational modulation of the Ho- 
line profile by the acc retion flow has pre viously been reported 
in another young BD (Scholz et al. 2005). A hot accretion spot 
would result in the periodic modulation of the light curve on 
a timescale of the BD rotational period. The 6h rise time of 
2MASS J0414's {/-band light curve suggests a rotational period 
of the order of 12 h, and is thus consistent with rotational mod- 
ulation. Assuming a black-body temperature of 8000 K for the 
accretion shock, the ho t spot would cover about 0.3% of the 
substellar surface ("e.g.. iBouvier et af This value is sim- 

ilar to the size of accretion shocks on classical T Tauri stars. 
The two subsequent U -band ground observations, obtained 6.4 h 
apart, show no sign of variability, unlike the OM light curve on 
this timescale. This indicates that 2MASSJ0414 was back to 
its quiescent level nearly one year after the OM observations. 
Accretion hot spots have been reported to be variable feat ures in 
T Tauri stars, and even in one young BD (see lScholz et alJ2005l 
and references therein). 

We conclude that at least two plausible sources of vari- 
ability could account for the observed f/-band light curve of 
2MASS J0414: an episodic accretion event onto the central ob- 
ject or the rotational modulation of the BD luminosity by a hot 
accretion spot. In both cases these observations indicate variable 
mass accretion rate in this BD. 

Our [/-band survey of 1 1 BDs (of which 5 are actively ac- 
creting and 4 are not accreting) indicates 2MASS J0414 as the 
brightest source, displaying variability related to accretion from 
a circumsubstellar disk. 2MASS J0414 is then the best target of 
this sample of 5 accreting BDs for ground follow-up observa- 
tions to investigate accretion in substellar objects, for example 
with simultaneous [/-band and spectroscopic ground monitoring 
with a one-hour sampling. 



6 The ballistic infall velocity from rest at t he disk's truncation ra - 
dius, r m , is: V(r) 2 = 2GM+ (l/r - l/r m ) fe.g.. lHartmann et aljfl994h . 
Therefore, assuming a purely radial infall, the travel time from the 
disk's truncation radius to the substellar surface is: t = r„ /V2GM* X 
(sin 26 -26 + n)/2, where 6 = arccos (VI - R-„/r m ). 
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Appendix A: Formulae for the OM aperture 
photometry and associated errors 

We give here the formulae for the canonical 6"-radius aperture 
that we used in our IDL procedures for the OM aperture photom- 
etry and associated errors. 

The relation between the count rate of incident photons, Co, 
and the measure count rate on the detector, Cdet, is: 

Co = F c f°(*, T) X G c e ; nP (x) X /deadnme^) X / tds , (A. 1) 

where F 'f 160 is the theoretical formula for coincidence-loss, G e , mp 

cl ' cl 

is the empirical correction of F* eo from calibration observa- 
tions, /deadtime is the deadtime correction, and fys is the temporal 
degradation sensitivity correction. T is the frametime, i.e., the 
sum of the frame transfer time (Tr = 0.1740ms), needed to 
transfer the charges accumulated on the continuously exposed 
CCD during an integration time into the storage area, and the 
time needed to read out all the pixels in the storage area. The 
frametime depends on the size and position of the exposure win- 
dow; here it is in units of seconds, and x = Cdet x T. Note that 
using an aperture radius different than the 6" value used in the 
calibration of the coincidence-loss would add a correction factor 
in Ea. lA.ll using the PSF calibration. 

The theoret ical correction of the coincidence-loss (see 
lEhle et al.l2005l) is given by: 

F$ eo {x,T)='^-£. (A.2) 
l ft - 1 

The empirical correction of this theoretical formula is: 

G™ p (x) = l+Z^fl,-*', (A.3) 

where the polynomial coefficients are provided by the 
extension COINCIDENCE of the current calibration file 
OM.PHOTTONAT_0003 .CCF. 

The only time when source photons are not properly 
recorded by the detector is during charge transfer; therefore the 
deadtime is equal to the frame transfer time, and the deadtime 
fraction is d = T^/T . Consequently, the deadtime correction is: 

/deadtime(r) = 1/(1 - d) = T/(T - 7ft). 

The temporal degradation sensitivity correction, / t d s , is avail- 
able in SAS version 6.5 in the fits header of the observation 
source list (*OM*SWSRLP). 

The count rate error of incident photons, ACo, is computed 
by Gaussian propagation in Eq. IA.ll of the Poissonian count rate 
error on the detector, ACdet = -y/Cdet/Texp, where r exp is the ex- 
posure time, which leads to: 

ACp = /tds T 2 
AC det (T - r ft ) 2 

xl j-^ ln(l-*)Zii*<»/* [• (A.4) 

We calculate the error associated with an object count rate 
measured with aperture photometry as the quadratic sum of the 
count rate error of incident sky+object photons and the count 
rate error of incident sky photons determined with previous 
equations. 
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